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Abstract 

The entropic corrections to the flux-line energy of extreme type-II superconductors are 
computed using a schematic dual Villain model description of the flux quanta. We find that 
the temperature profile of the lower-critical field vanishes polynomially at the transition 
with an exponent v = 2/3 in the isotropic case, while it exhibits an inflection point for the 
case of weakly coupled layers in parallel magnetic field. It is argued that vestiges of these 
effects have already been observed in high-temperature superconductors. 
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It is well known that type-II superconductors allow magnetic fields to penetrate in the 
form of flux lines quantized in units of the flux quantum, $ = hc/2e, for field values above 
the lower critical field given by H cl {T) = 47T£(T)/$ 05 where e(T) denotes the flux-line free 
energy per unit length. 1 The London approximation yields Sq(T) = (<&o / 4ir \l) 2 \ii \l / £ for 
this free energy, where both the London penetration length, A^, and the coherence length, 
£, vary as (T c0 — T) -1 / 2 near the critical temperature, T c0 , within the Ginzburg-Landau 
mean-field approximation. Spatial fluctuations of a flux line act to reduce this mean- 
field line tension, however, particularly near the transition temperature. Such entropic 
corrections to the lower critical field have been computed 2 ' 3 and observed 4-6 in the case of 
high-temperature superconductors, where fluctuation effects are important due to the short 
coherence lengths and the large transition temperatures that characterize these materials. 

In this paper, we shall re-compute the flux-line tension of extreme type-II supercon- 
ductors (Xl 3> in the absence of magnetic field utilizing the so-called "frozen" limit of 
lattice superconductors first studied by Peskin. 7-9 The principle idea here is to exploit the 
duality of quantized magnetic flux lines in superconductors to confining strings of electric 
flux in pure lattice gauge theories. Hence, only the entropic correction to the line tension 
due to flux-line wandering is accounted for. The energy functional for three-dimensional 
(3D) superconductors may be expressed as 7 

E = Ps £{1 - cos [ A m<KO - A,(r))} + i A »Mr) - A,A M (r)] 2 , (1) 

where 4>{r) is the phase of the order-parameter on a cubic lattice with spacing a' > £, 
and where the magnetic flux threading the plaquette at site r perpendicular to the [i = 
x,y,z direction reads ^(r) = ($o/27r) J" Also, A M 0(r) = 0(r + p.) - 

<fi(r) is the lattice difference operator. A "frozen" superconductor (FZS) is simply the 
thermodynamic state in the limit of infinite phase rigidity, p s — > oo. This means that the 
London penetration length is small in comparison to the effective lattice constant, a', and 
hence that flux lines do not interact. In particular, upon making the gauge transformation 
= A'^ + A^cj), we immediately see that this limit implies A'^ = 2nm^, where m^(r) is an 
integer field. Therefore, the energy functional (1) reduces to E FZ = J2 r M (2e ) -2 [27rn M (r)] 2 , 
where n M (r) is the integer field describing unit line segments of quantized flux traversing 
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the plaquette (r, n) and satisfying the continuity equation 




(2) 



i.e., n M = Yliv 7 e MI / 7 A I/ m 7 . The universality class of the transition into the normal state 
thus falls within that of the inverted 3D XY model. 7-9 Employing this schematic descrip- 
tion of flux lines in the absence of external magnetic field, with an effective lattice constant 
a' ~ \l, we shall demonstrate below that the lower critical field vanishes as H c i oc (T c —T) 1/ 
near the transition in "frozen" lattice superconductors, where the exponent v = 2/3 cor- 
responds to that of the 3D XY model. Via an anisotropic version of this FZS model, 10 
we shall also show that the temperature profile for the parallel lower critical field in super- 
conducting films with weakly coupled layers has an inflection point. Last, we argue that 
these effects are realized in high-temperature superconductors in the critical regime. 4-6 

Isotropic Superconductor. Let us then describe the zero-field thermodynamics of flux 
lines in 3D superconductors via the energy functional 



along with constraint (2), where n M (r) is the integer field that represents a unit line segment 
of magnetic flux quanta threading the plaquette (r, fi) . The energy scale above is set by the 
mean-field flux-line free energy, £o(^) = ($o/47rAo) 2 (l — T/T c0 ), and the effective lattice 
constant, a', yielding [3 = a'e (T)/kBT. Here, A is on the order of the zero-temperature 
London penetration length, while T c q denotes the mean-field critical temperature. We 
presume that the lattice constant satisfies a' ~ Al, which insures that the magnetic flux 
is quantized on the scale of an elementary plaquette. To compute the dimensionless flux- 
line tension, a, we shall place a unit monopole/anti-monopole pair separated by a large 
distance R ^> a' along a principle axis of the cubic lattice. The partition function for 
the FZS system (3) in the presence of such an external monopole charge configuration, 
p{r) = 5 r fi — 5 r} Rxi therefore reads 



Yet since Z[p]/Z[0] oc e ^ aR coincides with the correlation function of the dual 3D XY 




(3) 





model in the Villain form, 7 C(R) oc e the dimensionless flux-line free energy is then 



3 



related to the correlation length, £, of the latter by a = (Henceforth, R and £ 

are given in units of a'.) This implies that the line tension varies as a = cr (/3 — f3 c ) v 
near the critical temperature, where v = 2/3 is the critical exponent corresponding to 
the divergent correlation length of the 3D Villain model. 7 Here, (3 C = | is the critical 
temperature of the latter, 9 while cro is another non-universal constant of order unity. Since 
the dimensions of a typical flux loop are set by the divergent correlation length, £, short- 
range interactions among the flux lines on the scale of the finite mean-field penetration 
length should be irrelevant at criticality. This validates the FZS limit (3). The flux-line 
tension (or Z[p]/Z[0]) can also be computed at low temperature by permitting maximum 
transverse excursions in the line of one lattice spacing in any of the remaining 2(d — 1) 
directions, where d = 3 gives the dimension of the cubic lattice. This yields 

a = l-2(d-l)p- 1 e- f3 , (5) 

which agrees with the low-temperature dependence of the interface tension in the two- 
dimensional (2D) Ising model for the case d = 2, 11 as it should. 

The mean-field flux-line free energy in the "frozen" superconductor model [(2) and 
(3)] is therefore renormalized down by spatial fluctuations to e(T) = e (T)a(/3) oc (f3 — (3 c ) 1 ' 
near the critical temperature, T c = (T~ Q +[3 C T^~ 1 )~ 1 , which is determined by the condition 
(3 = (3 C . Here, the energy scale /cbT = a'($o/47rAo) 2 sets the size of the critical regime, 
5T C = T c0 —T c ~ T c 2 /T . In Fig. 1, we plot the predicted profile of the lower critical field vs. 
temperature, which is determined by the flux-line tension in zero-field via H c i = 47re/$o 
in the case of extreme type-II superconductors. Here, we have interpolated between the 
low-temperature and critical behaviors discussed previously. Last, we mention that a A- 
transition is also expected in the specific heat vs. temperature. 7 ' 9 

Layered Superconductor (Film). Consider now an anisotropic FZS model made up of 
iV weakly coupled square lattices separated by the same effective lattice constant a' ~ Ao, 
with an energy functional given by 

N N-l 

E FZ /k B T = pJ2/Z w l )\ 2 +^PzZH z )i 2 ( 6 ) 

1=1 r 1=1 r 

along with constraint (2). Here, n = (n x , n y ), r and / respectively denote the in-plane and 
layer indices, while 7 » 1 denotes the anisotropy parameter. 1 The smallest energy scale 
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is set by the mean-field in-plane flux-line tension, 1 r )~ 1 e (T), and the lattice constant, 
yielding (3 = / y~ 1 e (T)a / /ksT. The author has recently studied this model in the context 
of the dual XY model in the Villain form with a finite number of weakly coupled layers, 10 
where he finds a superfluid ordering transition at (inverted) temperature (3 C ~ 7r/4, followed 
by a layer-decoupling transition at higher (inverted) temperature /?* = 2tt. Very similar 
results were obtained previously by Korshunov in his approximate study of the infinite- 
layer case. 12 ' 13 Specifically, the in-plane field may be expressed as n(f, /) = n'(r,l) — 
n_(r, /) + n_(r, / — 1), along with constraints V • ft' = p and V • n_ = n 2 , where p{r) 
again represents the monopole distribution and V = (A x ,A y ). Employing the potential 
representation n_ = — V<E> factorizes the partition function (4) into Z = Zcg^iLiZ^ g 
in the limit 7 — > 00, 10 where Z^ G is the partition function for the 2D discrete gaussian 
model 14 that describes the Z th layer isolated from its neighbors (n z = 0), while the inter- 
layer Coulomb gas factor reads 

f N 

Z CG \p}= exp<^-/3^^K(f,/-l)-n z (f,/) + 2p(f,/)]G( 2 )(f-f / )x 

{n z (r,l)} { 1=1 r,r ' 

N-l ^ 

x [*,(?', I - 1) - n z (f', /)] - 7 /3 E n l<?> Z ) ' ( ? ) 

1=1 r J 

with the boundary layers set to n z (f, 0) = = n z (r,N). Here, = —V -2 is the 

Greens function for the square lattice. Notice that flux lines passing vertically through 
all of the layers make no contribution to the energy functional above. 15 In the limit of 
extreme anisotropy 7 — > 00, we thus find that Z^ G and Zcg respectively exhibit an 
inverted Kosterlitz-Thouless (KT) transition at (3 = (3 C and a Coulomb gas 16 transition at 
= 0*~ 10&. 

To compute the parallel flux-line tension, we again place unit oppositely charged 
monopoles separated by a large distance R ^> a' along a principle axis of a given layer 
/. Since the corresponding intra-layer correlation function of the dual Villain model fac- 
torizes into C(R) = C\\(R)C±(R), where C\\(R) = Z^ G [p} / Z^ G [0} = (r /R)^exp(-R/^) 
and C±(R) = Zcg[p]/-^cg[0] = (R/ro) v± , we obtain that the dimensionless parallel line 
tension is given by a = (/3^||) _1 , where £y denotes the correlation length of the 2D Villain 
model 16 dual to the former discrete-gaussian model. 10 ' 14 Hence, this line tension vanishes as 
a = a expl-Ao/iP-fic) 1 / 2 ] just below the critical temperature, 14 T c = (T^ 1 + P^Tq 1 )' 1 , 



where o"o and A again are non-universal constants of order unity. A similar temperature 
dependence for the parallel lower-critical field in layered superconductors has been obtained 
via a fermion analogy for the Lawrence-Doniach model. 13 The low-temperature behavior of 
the parallel line tension, on the other hand, is determined by Eq. (5) for d = 2. Although 
the corresponding parallel flux-line free energy, £\\(T) = r y~ 1 eo(T)a((3), is renormalized 
down with respect to the mean-field result just like in the isotropic case, we also observe 
that it must have an inflection point as a function of temperature just below the transi- 
tion, where it vanishes exponentially. This is displayed by the parallel lower-critical field 
(H^) in Fig. 2, which interpolates between the low-temperature and critical behaviors. 
The exponential suppression of near T c indicates that critical flux-line wandering is 
extremely violent in two dimensions. 

Concerning the lower-critical field perpendicular to the layers, classical electro- 
dynamics dictates that it must vanish in the present film geometry. However, we may fix 
the ends of a vertical flux line at the boundary layers [n z (r, 0) = <^o = n z(j-> N)] and com- 
pute its mean-squared displacement uf = J2f^o r2 { n zif, l)n z (0, 1))/ J2r^o( Uz (^ n z(0, /)) 
at a given layer /. In the limit of extreme anisotropy, 7 — > 00, it becomes sufficient to 
consider the double-layer case, N = 2. Since the charge correlation function for the 2D 
Coulomb gas (7) satisfies 16 (n z (0)n z (r)) oc -r" 2 ^, we obtain u 2 = a' 2 (/3- 
Exactly the same result is also obtained for pancake vortices in the absence of Josephson 
coupling between layers 3 ' 17 if the in-plane London penetration length is chosen for 
the lattice constant a'. The latter choice insures flux quantization at all temperatures. 
It also indicates that an elementary segment of vertical flux line, n z = ±1, corresponds 
to a stack of pancakes. We may therefore interpret the lower transition temperature, 
T* = (T^q 1 + P^Tq 1 ) -1 , determined by [5 = /3* = 2rr as an evaporation point for vertical 
flux lines (see Fig. 2), since u 2 diverges for T > T*. If 7 is now presumed to be large 
but finite, however, we must then require that the perpendicular magnetic field exceed the 
decoupling cross-over field at T*, 18 B^r ~ $o/7 2 <^ 2 - This guarantees that the distance be- 
tween neighboring stacks of pancake vortices be within the Josephson penetration length, 
7<i. Here d measures the distance between adjacent Josephson-coupled planes in the un- 
derlying Lawrence-Doniach model for layered superconductivity, 1 which is essentially an 
anisotropic version of Eq. (1). The former condition excludes the possibility that a well- 
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defined Josephson vortex connects displaced stacks of pancake vortices in adjacent layers, 
an effect that is not accounted for by the present factorization of the layered FZS model (6) 
into parallel and perpendicular parts. Last, we note that compatibility with the present 
FZS model requires that the physical length scales satisfy ^> a' ~ Ao- 

The above results for the FZS model (6) are in apparent contradiction with those of 
the more general anisotropic lattice superconductor model (1) obtained by Korshunov, 12 
who finds that layers decouple at a temperature tenfold greater than the superfluid KT 
transition. This paradox is resolved by observing (i) that the lattice constant in such 
models is on the order of the coherence length, £ <C Al, and (ii) that such decoupling of 
layers naturally coincides with the vanishing of the parallel lower-critical field here. Notice 
this means that in the regime T* < T < T c , a Meissner phase exists for in-plane magnetic 
fields while each layer is resistive, a feature which has been observed experimentally. 6 ' 13 
The latter, however, is then consistent with the "evaporation" of vertical flux lines at T* 
obtained here. Last, we mention that the line tension for parallel vortices can be computed 
directly from the double-layer XY model, 10 where one recovers the previous exponential 
temperature dependence shown in Fig. 2. 19 Again, this demonstrates the consistency of the 
present results with those obtained from model (1) in the extreme type-II limit, A^ — > oo, 
which is precisely the Xy-model. 

High- Temperature Superconductors. Let us now apply the above theory to the ox- 
ide superconductors, which are extreme type-II (Xl/C ~ 100) with observable fluctuation 
effects. 1-6 Consider first bulk YBa2Cu307_ y , which may be classified as a 3D superconduc- 
tor with strong anisotropy, 1 7 ~ 10. The scaled isotropic FZS model (3) predicts a width in 
temperature for critical fluctuations (see Fig. 1) on the order of 5T C ~ 7T C 2 /T , where again 
ksTo = ($ /47rAo) 2 a / . Assuming typical parameters T c0 ~ 100K and a' ~ A ~ 10 3 A, we 
obtain 5T C ~ 0.5 K. This estimate is consistent with recent H c i measurements performed 
by Pastoriza et al., 4 who observe fluctuation effects a la Fig. 1 to within this scale. If 
we suppose, however, that a' ~ A^(T) = Ao(T c q/5T c ) 1 / 2 , which insures flux quantization 
at all temperatures, then we obtain a much more modest estimate of 5T C ~ 2.5 mK for 
the size of the critical region. Also, the exponent v = 2/3 predicted here for the critical 
temperature dependence of the lower-critical field is consistent with the critical tempera- 
ture dependence of the London penetration length Al per the (neutral) 3D XY model; 1 
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i.e., H c i oc A^ 2 , with A^ oc (T c — T) -1 / 3 . The latter critical temperature dependence has 
been observed in recent penetration depth measurements on this class of high-temperature 
superconductors. 20 

The oxide superconductor Bi 2 Sr 2 CaCu208, on the other hand, may be classified as a 
layered superconductor with extreme anisotropy, 1 7 ~ 100. The layered FZS model (6) also 
predicts a fluctuation regime on the order of 5T C ~ 7T C 2 /T wide in temperature. Taking 
the previous values for A and T c0 , then the choice a' ~ A results in 8T C ~ 5 K, which is on 
the order of the effect (see Figs. 1 and 2) observed by Brawner et al. 5 and by Wan et al. in 
their measurements of the field of first penetration, H p , along the (perpendicular) c-axis for 
this material. On the other hand, the choice a' ~ \ L (T) = A (T c0 /5T c ) 1 / 2 for the lattice 
constant with the same parameter values results in a smaller estimate 5T C ~ 0.25 K for the 
size of the critical region. More interestingly, Wan et al. also observe an inflection point 
in Hp vs. T on this scale. This may be a vestige of the same prediction made here for the 
parallel lower-critical field if geometrical demagnetization effects are presumed to be strong. 
In particular, only parallel flux lines contribute to the total line free energy at temperatures 
T* < T < T c , since the perpendicular flux lines have by then "evaporated". This scenario 
requires, of course, that the perpendicular magnetic field exceed the decoupling cross- 
over scale, 18 B^r ~ $ /A 2 , where Aj represents the Josephson penetration length. Since 
this scale is known to diverge exponentially as T approaches the decoupling transition 
temperature T c from below, 19 the former cross-over scale should be negligibly small in 
that regime. At the evaporation point (T = T*), however, we have that Aj ~ 7<i. This 
implies a decoupling cross-over field of approximately 1 kG in the present case, 1 where 
d = 15 A. A direct measurement of the parallel lower-critical field in such material would 
clearly be very useful. Last, we mention that Wan et al. have fit their results for the 
2D-3D crossover line of c-axis aligned flux lines to 6 H m oc {(3 — f3 c ) u , with v = 0.68. The 
latter exponent is precisely that predicted for H cl vs. T in the isotropic FZS model (3), 
which suggests that H m is related in some way to bulk H cl along the c-axis. 

Finally, we note that the Coulomb gas model (7) for vertical flux-line segments in- 
dicates that the evaporation transition at f3~ x = f3~ x ~ / 9~ 1 /10 occurring in the limit 
7 — > 00 is renormalized down to f3~ x = e ^Z^ 1 f° r 7 finite, where eo — 1 oc e _27/? * 
is the polarization at the transition point due to flux-line/anti-flux-line excitations in be- 
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tween consecutive layers. Hence, flux line "evaporation" in clean layered superconductors is 
driven to lower temperature as the c-axis thickness increases. 21 If we presume the existence 
of strong columnar pins, however, then such isolated flux-line/anti-flux-line excitations be- 
tween layers are suppressed, resulting in an unrenormalized flux-line evaporation transi- 
tion. Precisely this effect has been observed in high-temperature superconductors, where 
the melting-evaporation-irreversibility line moves up in temperature upon the introduction 
of columnar pins. 22 

In summary, we have elucidated the critical behavior of the flux-line tension in both 
isotropic and layered extreme type-II superconductors (£ — > 0) utilizing the so-called 
"frozen" limit of a lattice model. 7-9 Notably, we predict an inflection point near the transi- 
tion for the parallel lower-critical field as a function of temperature in layered superconduct- 
ing films. Vestiges of this effect have already appeared in quasi-2D oxide superconductors. 6 
It is important to remark that fluctuations in the magnitude of the superconducting order- 
parameter are entirely neglected within this approach, contrary to a recent treatment of 
the same problem. 23 This is generally a valid approximation when effects due to the normal 
cores are negligible, like in the present limit £ — > 0, and in quasi-2D systems, where phase 
fluctuations drive critical properties. 

The author is indebted to H. Safar, M. Maley, S. Fleshier, G. Gilmer and L. Bu- 
laevskii for valuable discussions. This work was performed under the auspices of the U.S. 
Department of Energy, and was supported in part by National Science Foundation grant 
DMR-9322427. 
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Figure Captions 



Fig. 1. Shown is a temperature profile for H c \ (solid line) characteristic of isotropic extreme 
type-II superconductors in the critical regime (shaded area). Notably, H c i deviates 
from the mean-field result (dashed line) in an activated manner, and it approaches 
the transition temperature by following a polynomial law with exponent v = 2/3 (see 
ref. 7). 

Fig. 2. The parallel lower-critical field characteristic of extreme type-II superconducting films 
that are composed of weakly coupled layers is plotted as a function of temperature near 
the phase transition. Notice that it vanishes exponentially beyond the inflection point 
just below T c . Notice also that the low-temperature transition at T*, which marks 
where perpendicular flux lines "evaporate" , lies well inside the mean-field regime. 



12 



